Material dependence of Casimir forces: gradient expansion beyond proximity 
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A widely used method for estimating Casimir interactions [H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 
793 (1948)] between gently curved material surfaces at short distances is the proximity force approximation 
(PFA). While this approximation is asymptotically exact at vanishing separations, quantifying corrections 
to PFA has been notoriously difficult. Here we use a derivative expansion to compute the leading curvature 
correction to PFA for metals (gold) at room temperature. We derive an explicit expression for the amplitude 
Oi of the PFA correction to the force gradient for axially symmetric surfaces. In the non-retarded limit, the 
corrections to the Casimir free energy are found to scale logarithmically with distance. For gold, 6i has an 
unusually large temperature dependence. 



The continual drive towards miniaturization of de- 
vices inevitably leads to scales were quantum effects are 
significant. In the realm of micromechanical systems 
(MEMS)- nano fabricated devices actuated by electrical 
bias- Casimir^ and van der Waals forces^ are paramount. 
These forces originate in quantum fluctuations of elec- 
tromagnetic fields, and may cause devices to fail due to 
stiction."^ Due to their non-trivial dependence on shape 
and material properties, Casimir forces are difficult to 
compute.'*'^ A common method for estimating forces in 
a general setup is the proximity force approximation 
(PFA): First the force is calculated between two paral- 
lel plates at separation d with the Lifshitz formula (as 
a function of dielectric response of the adjoining slabs;^) 
then the geometry is taken into account by averaging 
over (appropriately defined) separations d between ad- 
joining surfaces.^ While this approximation is qualita- 
tively wrong in special circumstances,^ it remains a use- 
ful tool in high precision experiments between surfaces of 
large radii of curvature R. In these cases PFA is asymp- 
totically exact at small separations (for d <ti R). How- 
ever, improvements in sensitivity of measurements war- 
rant quantifying corrections to PFA which we undertake 
in this paper. 

We consider a geometry consisting of two infinitely 
thick plates composed of homogeneous and isotropic di- 
electric materials, with permittivities ei(w) and e2(w). 
For simplicity we assume that one of the plates is a plane- 
parallel slab, while the other is gently curved, and char- 
acterized by a smooth height profile z = -ff(x), where 
z is the local distance from the planar surface S2, and 
x = {xi,X2) is the vector spanning E2 (see Fig. 1).^ Fol- 
lowing Refs. 10 and 11 we postulate that the Casimir free 
energy J^[H] admits a local expansion of the form 

= J■pFA[i^]+ / d^a{H)VH-VH + p^'^\H], (1) 
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FIG. 1. Parametrization of a the profile of a gently curved 
dielectric surface near a flat dielectric plate. 

where Tpfa[H] represents the PFA free energy 

^pfa[H] - / d^TppiH) , (2) 

with J-pp{z) the free energy per unit area for two plane- 
parallel dielectric plates of permittivities ei(w) and e2('^) 
at distance z, as given by the Lifshitz formula,^ and a{H) 
is a function to be determined. The quantity p^^-* [H] in 
Eq. (1) represents corrections that become negligible if 
the curvature of the surface is small compared to the 
minimal distance d between the surfaces. "'^^ 

The function a{H) in Eq. (1) can be determined from 
a perturbative expansion of the Casimir free energy in 
the deformation profile /i(x), defined above the point of 
closest proximity as H{x) = d -\- h{x), see Fig. 1. Note 
that the latter perturbation requires a small deformation 
amplitude while the gradient expansion relies upon small 
changes in the slope. However, it can be shown that un- 
der certain conditions (existence of perturbation theory 
in ft.(x) and regularity of the involved kernels in momen- 
tum space) the gradient expansion of Eq. (1) follows from 
resumming the perturbative series for small in-plane mo- 
menta. If so, the expansion of the free energy in Eq. (1) 
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has to match the pcrtiirbativc expansion for |/i(x)|/d <C 1 
at small momenta. The latter is given by 



■/i(x)] = AF„(d) + ^l{d)h{Q 



G(fc;d)|Mk)p+p(2) 



[h], (3) 



where A is the surface area, k is the in-plane wave- vector, 
/i(k) is the Fourier transform of /i(x), and /0^^'[/i] refers 
to higher order corrections. The function a{H) can now 
be determined if the kernel G(k; d) can be expanded to 
second order in k (which as we shall see is the case for 
our problem). Indeed, matching the expansion 

G{k; d) = 7(rf) + 5{d) fc^ + o{k^) , (4) 

to Eq. (1) leads to 

T'^pid) = f,{d) , = 27(d) , a{d) = S{d), (5) 

where a prime denotes a derivative with respect to d. For 
further justification of the gradient expansion we refer to 
the discussion after Eq. (12). 

We now give an outline of the computation of the ker- 
nel G{k;d). The staring point is the scattering formula 
for the Casimir free energy J^,^^ 



T = )!cBr^'Trlog[l - T(i)UT(2)u] , 



(6) 



n>0 



where feg is Boltzmann's constant, T the temperature, 
and the primed sum runs over Matsubara frequencies 
^„ = I'KnksT jh with the n = term weighted by 1/2. 
In Eq. (6), T^^^ denotes the T-operator of plate j, evalu- 
ated for imaginary frequency In a plane- wave ba- 
sis |k, Q),^* where Q = E,M \^ the polarization, E 
and M denote respectively electric (transverse magnetic) 
and magnetic (transverse electric) modes. The transla- 
tion operator U in Eq. (6) is diagonal with matrix ele- 
ments e"'''" where g„ = ^k'^ + k\ = qn{k), k = |k| and 
Kn = £,n/c. For the undeformed slab, the operator T^^^ 
is diagonal in the plane-wave basis with matrix elements 
given by the scattering amplitudes 

rig,(k,k') = (27r)25(2)(k-k')<5QQ.rg)(ie„,k) , (7) 
where rQ^(i^„,k) are the Presnel reflection coeflicients 
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k"^ . There are no analytical for- 



and Sn^ 

mulae for the elements of the T-operator of the curved 
plate T(^), but for small deformations they can be ex- 
panded in powers of h{x) as^^ 



(k, k') = (2^)2^(2) (k - k') Sqq' rg' {zU, k) 



fd^k" 



-2i?QQ,(k,k')Mk-k') 
(B2)QQ-(k,k'; k")h{k-k")h{k" -k') + ... 



(9) 



where q'^ ~ qn{k'). The coefficients i3QQ'(k, k') and 
(i32)QQ(k', k'; k") can be obtained by standard pertur- 
bation theory in the height field and arc given in Ref. 15; 
they depend on the relative orientation of the wave vec- 
tors k and k', on the corresponding g„ and s„ and the 
dielectric function itself. Substituting the above expan- 
sion into Eq. (6), we obtain 

G(M) = tsTY,J§^, /„(k-.k-+k)+/„(k-.l.--l.) _ 

(10) 

where 

/„(k', k") = - X: HBMk', k'; k") 

+ 2E 4^^e-«:.''^Sgg,(k',k")Sg,Q(k",kO] , 

(11) 

with q'J^ = qn{k"), gQ{k) = 1 - rgVg^ exp[-2g„d]. The 
explicit dependence of several quantities on is not 
shown for brevity. Since the sum over Matsubara fre- 
quencies in Eq. (10) is exponentially convergent, the ex- 
istence of the second fc-derivative of G{k] d) is ensured 
if for all n > the second derivative of /„(k',k' -I- k) 
with respect to say kx, is absolutely integrable over k' 
for k = 0. We then have 
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Having obtained the amplitude function a{H) by a 
comparison with the pcrturbative free energy at small 
in-plane momenta, we note that in principle, the free en- 
ergy of Eq. (1) should follow directly from Eq. (6) when 
the T-operators are expanded in powers of the gradient of 
the surface profile /i(x), making no assumption about the 
amplitude of /i(x) itself. Indeed, such an expansion has 
been carried out for dielectric surfaces in Ref. 15. Substi- 
tution of this gradient expansion in Eq. (6) yields a non- 
local functional for the free energy. Hence, the function 
a{H) must be determined by the gradient expansion of 
the T-operator and subsequent locality expansion of the 
free energy. -'^^ 

The above formulae enable evaluation of a{d) for ar- 
bitrary dielectric functions ei{oj) and £2(0;). In the ideal 
limit of perfect conductors at zero temperature (as well as 
for scalar fields obeying Dirichlet, Neumann and mixed 
boundary conditions), a{d) has a simple power law de- 
pendence on d; as a{d) ~ 1/d^ with a coeflicient that 
can be computed exactly. This permits to obtain sim- 
ple closed formulae for the leading correction to the PFA 
for a variety of profiles h{x)}^ In general, for finite T 
and/or for dielectric materials, the dependence of a(d) 
on d cannot be expressed as a simple power law, and has 
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to be computed numerically. Evaluation of Eq. (1) then 
permits to estimate the leading correction to PEA for ar- 
bitrary shapes of the profiles, and for any materials and 
temperatures of the involved bodies. 

While numerical evaluation of Eq. (1) for specific ex- 
perimental setups (e.g. spheres, cylinders or corrugated 
plates) is straightforward in general, by simple manip- 
ulations of Eq. (1) not presented here for brevity it 
is possible to obtain a closed semi-analytical expression 
for the leading correction to PEA for the gradient of the 
Casimir force dF/dd = —d'^J^/dd'^ (a quantity that is 
measured directly in some experiments^^), in the exper- 
imentally relevant geometry of an arbitrary but axially 
symmetric profile /i(x), which is a smooth function of 

= xl+X2- We expand the profile around the point of 
minimum separation as 

2 4 

h{xi,X2) = h{p^) = d+^+Ci^ + ... , (13) 

where R sets the overall scale of curvature. By substi- 
tuting the above expansion into Eq. (1), one can derive 
the leading correction to PEA for the force-gradient as 

BF / d \ 

— = -27rRF„{d) [l + Oi^+ o{d/R)j , (14) 

where Fpp{d) = —dTpp{d) / dd is the Casimir force per 
unit area between two parallel plates, and 

Note that the dimensionless function /3(c?) depends on 
the ratio of d to one (or more) material dependent length 
scale. Independently of the dielectric material, the func- 
tion 5{d) becomes a simple power law in the non-retarded 
limit where J-pp{d) ~ 1/d^ and 5(d) ^ l/d^ so that f3{d) 
and hence 9i tends to a constant for d ^ 0. This implies 
that in the non-retarded limit the first correction to the 
PEA for J" scales logarithmically with the separation d. 

Using the above formulae, we numerically computed 
the coefficient 9i in the experimentally relevant sphere- 
plate geometry (corresponding to Ci = 1/4), assuming 
for simplicity that both surfaces are made of the same 
material; gold at T = 300 K. Eor a perfect reflector at 
zero temperature, an exact value for 6i can be obtained. 
In this limit the coefficient (3 was computed in Ref. 11, 
and given as /3em = 2/3(1 — IS/tt^). Since for perfect 
reflectors at T = 0, Tpp = -ir'^hc/ {12Qd^), it follows 
from Eq. (15) that ei\1f^ = i (2/3em - 1) = -0.564. 
Away from this ideal limit, 9\ has to be computed nu- 
merically. In our computations for gold, we adopted the 
Drude model (for a discussion of alternative models see 
Ref. 18, Chap. 13). In this model, the permittivity 
e(i^) to be used in the computation of the Casimir force 
for (ohmic) conductors is obtained using the well known 
Kramers-Kronig relation. For Im[e(a))] we used the tab- 
ulated optical data for gold quoted in Ref. 19, extrap- 
olated to zero frequency by means of the Drude model. 



with parameters fip = 9 cV /h and 7 = 35 meV/fi (see 
Ref. 18, p. 336). The resulting Ox is plotted in Eig. 2, as 
a function of logio(27rti/A3ooK), where At = hc/{kBT) 
is the thermal wavelength (A300K = 7.6 /im). The inset 
of Fig. 2 shows these data versus the separation in lin- 
ear scale (in microns). The crosses in Eig. 2 correspond 
to data for T — 300 K, while the dashed lines are for 
T = K. We also show (solid lines) the result for per- 
fect reflectors at T = 300 K. Thus, the solid lines neglect 
finite conductivity corrections, while the dashed lines ne- 
glect finite temperature corrections. The importance of 
including both corrections for the leading (PEA) Casimir 
force has been previously noted^°. A comparable previ- 
ous estimate of 9\ , for a metal sphere-plate setup at zero 
temperature, was obtained in Ref. 21. This estimate is 
based on a numerical summation of a partial wave series, 
valid at large separations d^ R, and fitting the results 
at short separations to a polynomial scries. While 
our explicit computation here confirms a linear correction 
to PEA as in Eq. (14), we expect higher order corrections 
to involve non-polynominal logarithmic terms^^. 

Some comments about the qualitative behavior of the 
numerical results are in order. As expected, the T = 
perfect conductor limit of -0.564 is approached by the 
perfect conductor data at T = 300 K (solid lines) for 
short separations, which is where temperature effects be- 
come negligible; while the T = data for gold (dashed 
lines) approach this limit for large separations, which is 
where the skin depth becomes negligible. We observe also 
that in the limit of large separations, the T = 300 K data 
in Eig. 2 (crosses and the solid line), both approach the 
limit 9i{d -J> 00) = 1/(12 C(3)) = 0.0693, which coincides 
with the classical high-temperature limit of a scalar field 
obeying Dirichlet boundary conditions. This is plausi- 
ble since for finite T in the limit of large separations the 
Casimir force is dominated by classical thermal fiuctua- 
tions of the electromagnetic field. It is well known that 
for non-magnetic materials like ohmic conductors, trans- 
verse electromagnetic fields decouple in the classical limit 
(Bohr- van Leeuwen theorem^^), and the Casimir force 
arises entirely from low-frequency fluctuations of the lon- 
gitudinal electric held (non-retarded limit of E modes), 
associated with scalar potentials which vanish on the sur- 
face of a conductor. 

Importantly, we see from Eig. 2 that the data for gold 
at T = 300 K deviate significantly from both the solid 

and the dashed curves for separations d between 150nm 
and 300nm, showing that in this range a precise de- 
termination of 9i requires simultaneous consideration of 
both finite conductivity and finite T corrections. This 
is somewhat surprising since the thermal correction to 
the Casimir force between two metallic plates is small for 
such small separations (in the Drude model, less than 3% 
below 200nm). On the contrary, we see from Eig. 2 that 
the magnitude of the thermal correction to 9i is about 
20% for d = 200nm, indicating that 9i is much more sen- 
sitive to temperature than the Casimir force itself. Eor a 
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FIG. 2. 9i for a gold sphere in front of a gold plate, computed 
using the tabulated optical data for gold, see Ref. 19. Crosses 
correspond to T = 300 K, while the dashed line is for T — OK. 
The solid line is for an ideal conductor at T = 300 K. For 
ideal conductors at T = 0, = —0.564 independently of 
separation. The inset depicts the same data, for gold at 300 K, 
as a function of separation (in microns) on a linear scale. 

vanishing separation we get 9i{d — > 0) = —0.206 which is 
consistent with our general prediction of a finite 9i in the 
non- retarded limit, and suggests that 6i(d) is a bounded 
function. 

Since Casimir force measurements are often limited to 
rather short separations between weakly curved surfaces, 
our work has immediate experimental relevance. An at- 
tempt to measure the coefficient 9i for a gold plate-sphere 
geometry at room temperature was made by using a mi- 
cromachined torsional oscillator. The best estimate of 
9i was obtained by dynamically measuring the Casimir 
force gradient at fixed separation d, for five spheres of 
different radii R, ranging from 10.5 to 148.2/im, and fit- 
ting the data linearly in 1/R. It was estimated that 9i is 
less than 0.4 at 95% confidence level, in the range 164nm 
< d < 300nm. Our data in Fig. 2 are slightly above, but 
roughly consistent with this bound. The small disagree- 
ment between our estimate of 9i and the experimental 
bound may be explained by observing that the theoret- 
ical prediction of 9i is sensitive to the plasma frequency 
ujp, with smaller ujp leading to smaller values for 9i. In 
our computations we used the standard value for gold 
{up = 9 eV/h), but it is well known that gold films pro- 
duced by deposition techniques may have much smaller 
values of the plasma frequency.^"' Another consideration 
is that, when fitting the force-gradient data in Ref. 23 
higher order corrections, of magnitude 92{d/R)^, to PFA 
were neglected. The coefficient 92 cannot be computed 
analytically yet, but its magnitude was estimated for 
perfect conductors at zero temperature through a Fade 



extrapolation constrained by a multipole expansion at 
large separations and by the gradient expansion at short 
separation. It appears that 02 is positive and of order 
one. Neglect of this second order correction in the fits of 
Ref. 23, especially for the smallest sphere, may lead to 
a systematic underestimation of 9i . Finally, it would be 
important to investigate the infiuence of surface rough- 
ness on 9i. 
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